Axisymmetric gravitational MHD equilibria in the presence of 

plasma rotation^ 

C. Cremaschini" , A. Beklemishev'', J. Miller'^''^ and M. Tessarotto'^'-'^ 
"'Department of Astronomy, University of Trieste, Italy, 
^ Budker Institute of Nuclear Physics, Novosibirsk, Russia, 
^ . International School for Advanced Studies, SISSA and INFN, 

^ \ Trieste, Italy, Department of Physics (Astrophysics), 

C I University of Oxford, Oxford, U.K., 

^ ' '^Department of Mathematics and Informatics, University of Trieste, Italy, 

^ . f Consortium of Magneto- fiuid- dynamics. University of Trieste, Italy 

(Dated: June 30, 2008) 

O ■ Abstract 

in : 

c/3 , In this paper, extending the investigation developed in an earher paper (Cremaschini et al., 2008), 

we pose the problem of the kinetic description of gravitational Hall- MHD equilibria which may arise 
^ ■ in accretion disks (AD) plasmas close to compact objects. When intense EM and gravitational 

cn : 

CN , fields, generated by the central object, are present, a convenient approach can be achieved in the 

a^ ■ 

\ context of the Vlasov-Maxwell description. In this paper the investigation is focused primarily on 

\o '• 

■ the following two aspects: 

oo ; 

. 1) the formulation of the kinetic treatment of G-Hall-MHD equilibria. Based on the identification 

' of the relevant first integrals of motion, we show that an explicit representation can be given for the 

' equilibrium kinetic distribution function. For each species this is represented as a superposition of 

suitable generalized Maxwellian distributions; 

2) the determination of the constraints to be placed on the fluid fields for the existence of the 

kinetic equilibria. In particular, this permits a unique determination of the functional form of the 

species number densities and of the fluid partial pressures, in terms of suitably prescribed flux 

functions. 

PACS numbers: 52.25.Dg,95.30.Qd,97.10.Gz 
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I. INTRODUCTION 



An open issue in astrophysics is understanding the dynamics of accretion disk (AD) 
plasmas occurring near compact objects (black holes, neutron stars, etc.) 0, [s, ^. 

It is well known that this problem must be generally formulated in the context of kinetic 
theory for the proper description of macroscopic plasma dynamics. In fact, unless direct 
information about the fluid fields which characterize the plasma is available (for example 
from theory, particle numerical simulations or direct observations), a purely fluid treatment 
of these plasmas may pose serious difficulties. These are related to the notorious impossibility 
of uniquely defining consistent closure conditions for the fluid equations. A convenient 
solution can, however, be reached by adopting a kinetic approach. An important motivation 
for this investigation is, in fact, that of avoiding the closure problem characteristic of fluid 
equations (see the related discussion in Paper I, Cremaschini et al., 2008). A flrst step in 
this direction, realized in this paper, is related to the investigation of kinetic equilibria (in 
this context, what is meant by an "equilibrium" is stationary-flow solution). The purpose 
of this paper is to investigate, in particular, kinetic axi-symmetric gravitational Hall-MHD 
(G-Hall-MHD) equilibria occurring in AD plasmas arising around compact stars, which are 
locally characterized by the presence of a family of nested axi-symmetric toroidal magnetic 
surfaces {ipir) = const} , ip denoting the so-called poloidal flux. These plasmas are expected 
to be coUisionless and characterized by the presence of intense EM flelds, as well as - at the 
same time - by a strong gravitational fleld generated by the central object. Extending the 
investigation developed in an earher paper (see Paper I) and based on a fluid description, in 
this paper we intend to analyze, in particular, the kinetic constraints placed on the relevant 
fluid fields by the kinetic treatment of G-Hall-MHD equilibria, i.e., the requirement that the 
kinetic distribution function is itself a stationary solution of the relevant kinetic equation. 
Ignoring possible weakly-dissipative and time- dependent effects as well as assuming the 
possible presence of a radial flow velocity in the AD, we shall assume - in particular - that the 
kinetic distribution function and the EM flelds associated with the plasma obey the system 
of Vlasov-Maxwell equations. In addition, the effect of the gravitational fleld produced 
by the central object is treated by the introduction of a pseudo-Newtonian gravitational 
potential which permits to retain weakly relativistic effects on the AD plasma dynamics. In 
such a case, in general, the form of the equilibrium kinetic distribution remains completely 
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arbitrary. The only restriction on its form, besides its (strict) positivity and the assumption 
of its suitable smoothness in the relevant phase-space, is obviously due to the requirement 
that it must be a function only of the independent first integrals of motion. As such, 
it is always possible to represent it as superpositions of suitable generalized Maxwellian 
distributions (GMD). In this paper we wish to analyze the constraints placed on them by 
the assumption of the existence of kinetic equilibria of this type. We intend to prove, in 
particular, that: 

1 . The assumption of kinetic equilibrium determines uniquely the functional form of the 
number densities, flow velocities and temperatures carried by each GMD, which are 
found to depend on appropriate flux functions to be suitably prescribed. 

2. In particular, we intend to prove that the flow velocities carried by each GMD are 
species-dependent but coincide, for each GMD belonging to the same species. 

3. In addition, the related angular velocities [Qg] are constant on each toroidal magnetic 
surface {'^(r) = const}, while the form of fig is uniquely prescribed in terms of nor- 
mal derivatives of suitable flux-surface averages of the electrostatic and gravitational 
pseudo-Newtonian potentials. 

II. KINETIC G-HALL-MHD EQUILIBRIA 

As a starting point, we require that AD plasma admits fluid equilibrium, i.e., a G- 
Hall-MHD equilibrium in the sense defined in Paper I. In particular, we assume that the 
equilibrium magnetic field B admits, at least locally, a family of magnetic surfaces = 
{■0(r) = const} , all mutually nested and bounded, which are represented by smooth toroidal 
surfaces, ip denoting the poloidal magnetic flux of B (see Paper I). By assumption, these 
surfaces and all equilibrium fluid flelds are axi-symmetric. This means that the relevant 
dynamical variables characterizing the plasma arc required to be independent of the toroidal 
angle when referred to a set of cylindrical coordinates {R,ip,z). In particular, the axis 
R — can be identifled with the principal axis of the toroidal surfaces. The assumption 
of the existence of a family of nested magnetic surfaces {ip{r)} implies that locally a set of 
magnetic coordinates (f, •&) can be deflned, with denoting, in particular, a curvilinear 
angle-like coordinate on a magnetic surface ip{r) — const. For deflniteness we shall assume, 
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furthermore, that the plasma is: a) non-relativistic, in the sense that - for an inertial observer 
at rest with respect to the central object and in the considered region of the AD - the 
species flow velocities are much smaller, in magnitude, than the speed of light in vacuum; 
b) collisionless, so that the mean free path of the particles of the plasma is assumed to be 
much larger than the largest characteristic scale length of the plasma; c) the plasma can be 
treated in the pseudo-Newtonian approximation, whereby the gravitational force produced 
by the central object is treated by means of an appropriate pseudo-Newtonian gravitational 
potential acting on the plasma. 

Let us now introduce the key assumption that the G-Hall-MHD equilibrium also corre- 
sponds to a kinetic equilibrium, namely that the kinetic distribution functions characterizing 
the plasma species arc stationary solutions of the Vlasov kinetic equation. We stress that 
the condition of kinetic equilibrium is here assumed to apply perhaps only in a local sense, 
i.e., in an appropriate sub-domain of the AD where the plasma can be treated, in particular, 
according to the previous assumptions a)-c). In particular, for greater generality and in 
contrast with respect to Paper I, here we shall allow that the equilibrium magnetic field 
admits a non- vanishing toroidal component (By), possibly produced by currents (of the AD 
plasma) located far from the considered region. This means that the kinetic distribution 
functions must be an exact first integral of motion. For definiteness we shall consider here 
a plasma formed by at least two species of charged particles: one species of ions and one of 
electrons. Regarding the specific form of the species-equilibrium kinetic distribution func- 
tions fk (where k — i,e in the case of a two-species plasma formed by one species of ions 
and the electrons), we shall assume that they can be approximated by a superposition of 
suitable GMD's f^g, namely for A; = i, e 

s£lk 

where /j and le are suitable sets of indices. Thus, in this sense for the kinetic description 
of the plasma, several sub-species s — 1, ..,n can in principle be distinguished, each one by 
assumption defined in such a way that: 

I) the kinetic distribution /*s is a first integral of motion, in the sense that: 

^ln/*. = 0. (2) 

II) for each s — l,n the equilibrium kinetic distribution function, denoted by /*s, is - in 
a suitable asymptotic sense to be defined later - a drifted local MaxweUian distribution of 
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the form: 



n, I M,fv-V,^^ 



This means that there exists an infinitesimal dimensionless parameter e such that 

hs = fus [1 + 0(e)] . (4) 
Here the notation is standard, thus n^, and are respectively the sub-species number 

1 /2 

density, temperature and flow velocity. Moreover, vths = {ST^/M^} is the thermal velocity 
and b is the unit vector of the magnetic fleld B(r,t). 

Ill) each distribution function J'ms {ioi s = 1, ra) is an approximate stationary solution of 
the Vlasov equation, in the sense that, consistent with the requirement (jlj), it must give 

^ln/Ms = + o(£). (5) 

It can be proved, as a consequence, that the following kinetic constraints must be imposed 
on the fluid flelds: 

1. Vs(r) can be be identified with a toroidal fiow velocity 

V, = e^m,, (6) 

where the toroidal angular frequency takes the form Qg = fig (■?/'), i.e., it is constant on 
magnetic surfaces ipi'"^) = const; 



2. the number density is: 



where 



Ss = Zse<j) + UGrav-^V,^ (8) 

and A = A — (A) . Here the notation is standard (see Paper I). In particular, (p and 
Ucrav are respectively the equilibrium electrostatic potential and an effective Pseudo- 
Newtonian gravitational potential. Moreover, (A) is the ^z^— surface average (defined 
on a fiux surface ipi^) = const.) of function A{r), by assumption independent of the 
toroidal angle ip, which is defined as (A) = f di}A{r)/ |B ■ V??! , with C, denoting 
^ = fdi9/\B-V^\ ; 
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3. ignoring possible slow spatial dependencies, the pseudo-density nos and the tempera- 
ture Tg are flux functions, namely they depend only on the poloidal flux ip so that: 

nos = riosi'ip) , (9) 
Ts = r.(^); (10) 

4. the toroidal angular velocity flsi4') can always be identified with the flux function 

nsW = n{^) + An,{^), (ii) 

where i^{ip) = c-^ ($) and A^ls{ip) = {Ucrav) are the toroidal angular velocities 
driven, respectively, by the electrostatic and effective gravitational potentials. 

Let us first determine the equilibrium distribution function f^^g- For this purpose, we first 
notice that it must be necessarily a function only of the independent first integrals of motion 
and/or of the relevant independent adiabatic invariants (at least for the leading orders in 
e). In axisymmetry these are well known j?, 8] in particular the canonical momentum p^g 
and the total particle energy. Eg = + ^$ + jUgrav In the following we assume that 
the plasma is strongly magnetized in the sense that for all species s = e,i 

"f-e, (12) 

where r^g = v±ths/^cs is the species average Larmor radius and Qcs = ZgcB/MgC is the 
species Larmor frequency. Here we shall assume that the plasma is subject to intense 
electrostatic and gravitational fields (defined respectively by the electrostatic potential $ 
and the effective gravitational potential Ugrav), thus requiring 

Tk^ir^^ (13) 

As a consequence, the magnetic flux as well as $ and the effective gravitational potential 
Ugrav are all considered to be of order while the species thermal velocities, i.e., v\\ths and 
vuhs (for 5 = 6,-?), are assumed to be of order o(e°). Constructing the adiabatic invariants 
in the following we shall consider in particular "thermal" test-particles for which |v| is of 
the same order as Vths (for s = l,n). It follows that, thanks to axisymmetry, the canonical 
momentum conjugate to the ignorable (toroidal) angle ip is by assumption a first integral. 
Taking into account the previous asymptotic orderings, it is given by 

p^g = MgRw-e^ + = —'ip^g (14) 

ec ec 



Similarly at equilibrium the test particle total energy, Eg = + ^Ugrav is by 

assumption a first integral of motion. It follows that the dynamical variable 

is manifestly a first integral of motion too. We notice that introducing the velocity vector 

= e^RQs{ip^s), it follows that H^s is also given by H^^ = ERsiVs) - ^i^^s{i^*s) - 
where 

so that E[is(Vs) can be interpreted as the total energy relative to the reference frame mov- 
ing with velocity V^.^. Next, we note that invoking the orderings f[T^ . f[T^ and assuming 
|v| /vths ~ the variables il'*s and if*^ can be conveniently approximated as 

= ^p[l + o{e)] (16) 

= Hs[l + oie^)] (17) 



where 



and Episi^s) is the total energy relative to the reference frame moving with the (fluid) 
velocity = e^RQsW [see Eq.(l6l)]. 

Thus, a possible definition for the equilibrium distributions /^.^ is provided by the distri- 
bution function 

{generalized Maxwellian distribution), where n*^ and T^^ are defined as n*^ = ns{ip*s) and 
T^.s = Ts{ip^s)- It clearly follows that: 

• for all s = l,n, f^,s is a function of first integrals and hence, by construction, is itself 
a first integral of motion. This result holds, in principle, for an arbitrary function 
^s{i'*s) (assumed to be suitably smooth); 

• the kinetic constraints fp])- ffTT]) are necessarily identically satisfied, by suitable defini- 
tion of the functions ns{ip*s) and fls{4'*s)] 

• invoking the asymptotic orderings (fT2|) . (fT3|) and |v| /vths ~ o(e°), Eq. (fT7|) is satisfied 
for Jms defined by Eq.([3]). 



To prove these statements, we invoke Eqs. (fT6l) and (fTTI) so that /*s can be written 

By making use of the expressions 

nosiij) = n,(V^) exp (-^ j , (21) 



s 



it follows that Eq.([2U]) leads to (HD, while Eq.([5I]) implies ([7]). Finally, we notice that 
Eq. ffTTl) can also be obtained in an equivalent way from the momentum balance equation 
for the s— th species [see Eqs.(l) and (2) in Paper I)]. In the present case, this gives 



qsfh 

e 



(^-V^+E^'^)) - j-nsVs X B = 0, (22) 



where Es is a suitable diamagnetic electric field (see Eq.(3), Paper I). It follows = 

a consequence, taking the scalar products term by 

term with b, and Vip it follows that 

V\\ = QR^, (23) 

l]i?|| = V^-(V0-Ei^)) (24) 
(-V0+Ei'^)) -B, = , (25) 

which proves Eq. (fTT|) . This result shows that in such for kinetic G-Hall-MHD 

equilibria which hold under the validity of the previous asymptotic orderings, no bifurcation 
arises (in contrast to the general case of fluid equilibria, considered in Paper I). This is due to 
the functional form of the number density provided by ([7j). As a consequence, as indicated 
above, in the present case the fluid velocity (V^) is uniquely determined and is identical 
for each sub-species belonging, respectively, to the ion or electron species. In addition, the 
presence of a finite (equilibrium) toroidal magnetic field is in principle permitted. Although 
its possible origin has not been explicitly discussed here, we can always assume that it 
is produced by "external" currents (i.e., located outside the domain of local existence of 
the kinetic equilibrium). Finally, another interesting consequence concerns the form of the 
generalized Grad-Shafranov equation, which in the present case becomes 



1 d^i/j 1 d^ijj _ AttR 
R~dz^ ^ RdR^ ~ c~ 



(26) 
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This shows that in kinetic G-Hall-MHD equihbria, the electric current density, producing 
the self-generated poloidal magnetic field, contains two terms: a) the first one driven by the 
'j/'-derivative of the average electrostatic potential is proportional to the local charge density 
of the plasma p (Hall effect); b) the second one due to the gravitational potential, occurs 
even in the case of a quasi-neutral plasma (p = 0). 



III. CONCLUSIONS 



In this paper, kinetic G-Hall-MHD equilibria have been investigated. In particular, the 
kinetic equilibrium has been determined imposing three main assumptions: 

1) that a fluid G-Hall-MHD equilibrium exists which is characterized locally by a family 
of nested magnetic surfaces {V'(r) = const} represented by axi-symmetric nested tori; 

2) requiring that the species equilibrium kinetic distribution function can be represented 
by a superposition of suitable local equilibrium distributions /*s (for s = l,n); 

3) imposing that the equilibrium distribution function can be expressed in terms of flrst 
integrals of the motion. 

Regarding the second assumption we remark that: a) experimental observations of col- 
lisionless astrophysical plasmas (for example, the solar wind) are consistent with this type 
of representation; b) it is always possible to represent a kinetic distribution function by an 
appropriate superposition of drifted local Maxwellian distributions (here replaced by the 
equilibrium distributions f^s)- Nevertheless, as shown in this paper, the assumption of ki- 
netic equilibrium invoked here has been obtained by imposing the validity of suitable kinetic 
constraints on these distributions. As in Paper I, no assumption of local quasi-neutrality has 
been invoked, thus permitting the consistent treatment of Hall effects. In addition, allowance 
has been made for the inclusion of relativistic effects, taken into account by means of an 
effective gravitational potential (Ugrav), and the presence of a flnite toroidal magnetic fleld. 
Finally, we remark that the third assumption actually rules out the treatment of possible 
more general kinetic equilibria, which may be expressed in terms of adiabatic invariants to be 
established based on the so-called gyrokinetic theory of particle dynamics in strong EM flelds 
The main results of the present theory include: 1) the treatment of kinetic equilibria 
in the presence of strong EM and gravitational flelds, i.e., in particular within the validity 
of the asymptotic orderings (|T2l) . (|T3l) : 2) the evaluation of the species-dependent toroidal 



9 



angular velocity (Qi), which - to leading order in the asymptotic parameter e, is proved to 
depend on the gradient of the ?/^— average electrostatic and effective gravitational potentials; 
3) the proof that for this type of kinetic equilibrium, no bifurcation occurs in the ion toroidal 
angular velocity {^i). The theory developed here applies to plasmas subject to the action 
of an intense magnetic field, partly self-generated by the plasma itself. The theory permits 
also the treatment of intense electric and gravitational fields, including the possible action 
of relativistic effects produced by the presence of a massive central object. In particular, the 
equilibrium distribution f^g has been determined by adopting a gyrokinetic approach. This 
permits f^^s to be constructed in such a way that it is - at the same time -a first integral of the 
motion and is also close, in an appropriate asymptotic sense, to a drifted local Maxwellian 
distribution. Based on the present theory, G-Hall-MHD equilibria can be investigated, in 
principle, utilizing either standard numerical solution methods or by a perturbative solution 
method based on a power series expansion near to the equatorial plane {z = 0). Detailed 
results will be presented elsewhere [10|. The results presented in this paper appear to be 
potentially significant both for the numerical evaluation of AD plasma equilibria and for the 
physical interpretation of the theory. We believe, in fact, that applications of the present 
theory can be useful for the theoretical investigation of AD phenomenology and for gaining 
a better understanding of the related basic physical processes. 
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